/*
 * @Descripttion:
 * @version:
 * @Author: xiaozilai
 * @Date: 2022-11-03 10:09:52
 * @LastEditors: xiaozilai
 * @LastEditTime: 2022-11-03 16:23:47
 */
/*
 * @lc app=leetcode.cn id=1620 lang=cpp
 *
 * [1620] 网络信号最好的坐标
 *
 * https://leetcode.cn/problems/coordinate-with-maximum-network-quality/description/
 *
 * algorithms
 * Medium (37.50%)
 * Likes:    87
 * Dislikes: 0
 * Total Accepted:    26.3K
 * Total Submissions: 56.5K
 * Testcase Example:  '[[1,2,5],[2,1,7],[3,1,9]]\n2'
 *
 * 给你一个数组 towers 和一个整数 radius 。
 *
 * 数组  towers  中包含一些网络信号塔，其中 towers[i] = [xi, yi, qi] 表示第 i 个网络信号塔的坐标是 (xi, yi)
 * 且信号强度参数为 qi 。所有坐标都是在  X-Y 坐标系内的 整数 坐标。两个坐标之间的距离用 欧几里得距离 计算。
 *
 * 整数 radius 表示一个塔 能到达 的 最远距离 。如果一个坐标跟塔的距离在 radius
 * 以内，那么该塔的信号可以到达该坐标。在这个范围以外信号会很微弱，所以 radius 以外的距离该塔是 不能到达的 。
 *
 * 如果第 i 个塔能到达 (x, y) ，那么该塔在此处的信号为 ⌊qi / (1 + d)⌋ ，其中 d 是塔跟此坐标的距离。一个坐标的 信号强度
 * 是所有 能到达 该坐标的塔的信号强度之和。
 *
 * 请你返回数组 [cx, cy] ，表示 信号强度 最大的 整数 坐标点 (cx, cy) 。如果有多个坐标网络信号一样大，请你返回字典序最小的 非负
 * 坐标。
 *
 * 注意：
 *
 *
 * 坐标 (x1, y1) 字典序比另一个坐标 (x2, y2) 小，需满足以下条件之一：
 *
 *
 * 要么 x1 < x2 ，
 * 要么 x1 == x2 且 y1 < y2 。
 *
 *
 * ⌊val⌋ 表示小于等于 val 的最大整数（向下取整函数）。
 *
 *
 *
 *
 * 示例 1：
 *
 *
 * 输入：towers = [[1,2,5],[2,1,7],[3,1,9]], radius = 2
 * 输出：[2,1]
 * 解释：
 * 坐标 (2, 1) 信号强度之和为 13
 * - 塔 (2, 1) 强度参数为 7 ，在该点强度为 ⌊7 / (1 + sqrt(0)⌋ = ⌊7⌋ = 7
 * - 塔 (1, 2) 强度参数为 5 ，在该点强度为 ⌊5 / (1 + sqrt(2)⌋ = ⌊2.07⌋ = 2
 * - 塔 (3, 1) 强度参数为 9 ，在该点强度为 ⌊9 / (1 + sqrt(1)⌋ = ⌊4.5⌋ = 4
 * 没有别的坐标有更大的信号强度。
 *
 * 示例 2：
 *
 *
 * 输入：towers = [[23,11,21]], radius = 9
 * 输出：[23,11]
 * 解释：由于仅存在一座信号塔，所以塔的位置信号强度最大。
 *
 * 示例 3：
 *
 *
 * 输入：towers = [[1,2,13],[2,1,7],[0,1,9]], radius = 2
 * 输出：[1,2]
 * 解释：坐标 (1, 2) 的信号强度最大。
 *
 *
 *
 * 提示：
 *
 *
 * 1 <= towers.length <= 50
 * towers[i].length == 3
 * 0 <= xi, yi, qi <= 50
 * 1 <= radius <= 50
 *
 *
 */

// @lc code=start
#include <iostream>
#include <vector>
#include <cmath>
using namespace std;
class Solution {
public:
    vector<int> bestCoordinate(vector<vector<int>> &towers, int radius) {
        buildSearchRange(towers);
        vector<int> ans{0, 0};
        int max_power = 0;
        for (int x = search_min_x; x <= search_max_x; ++x) {
            for (int y = search_min_y; y <= search_max_y; ++y) {
                int single_power = culSignalPower(towers, radius, x, y);
                if (single_power > max_power) {
                    max_power = single_power;
                    ans[0] = x;
                    ans[1] = y;
                }
            }
        }
        return ans;
    }

private:
    int search_min_x;
    int search_min_y;
    int search_max_x;
    int search_max_y;
    void buildSearchRange(vector<vector<int>> &towers) {
        search_max_x = towers.front()[0];
        search_min_x = towers.front()[0];
        search_max_y = towers.front()[1];
        search_min_y = towers.front()[1];
        for (auto &tower : towers) {
            search_max_x = max(search_max_x, tower[0]);
            search_max_y = max(search_max_y, tower[1]);
            search_min_x = min(search_min_x, tower[0]);
            search_min_y = min(search_min_y, tower[1]);
        }
        return;
    }
    int culSignalPower(vector<vector<int>> &towers, int radius, int x_pos, int y_pos) {
        int ret = 0;
        for (auto &tower : towers) {
            int tower_x = tower[0];
            int tower_y = tower[1];
            int tower_power = tower[2];
            double dist = culDistance(tower_x, tower_y, x_pos, y_pos);
            if (dist <= radius) {
                ret += (int)floor(tower_power / (1 + dist));
            }
        }
        return ret;
    }
    double culDistance(int tower_x, int tower_y, int x_pos, int y_pos) {
        return sqrt(pow((tower_x - x_pos), 2) + pow((tower_y - y_pos), 2));
    }
};
// @lc code=end
